Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 September 2009 (MN M£X style file v2.2) 



Can magnetic fields be detected during the inspiral of 
binary neutron stars? 

Bruno Giacomazzo, 1 Luciano Rezzolla 1 ' 2 and Luca Baiotti 3 

1 Max- Planck- Institut fur Gravitationsphysik, Albert-Einstein-Institut, Potsdam- Golm, mi6, Germany 

2 Department of Physics and Astronomy, Louisiana State University, Baton Rouge, LA 70803, USA 

' 5 Graduate School of Arts and Sciences, University of Tokyo, Komaba, Meguro-ku, Tokyo, 153-8902, Japan 



o 
o 

(N 
Oh 



o 
cr 



(N 
> 

(N 

(N 

o 

On 
O 



X 



7 September 2009 



ABSTRACT 

Using accurate and fully gcncral-rclativistic simulations we assess the effect that mag- 
netic fields have on the gravitational-wave emission produced during the inspiral and 
merger of magnetized neutron stars. In particular, we show that magnetic fields have 
an impact after the merger, because they are amplified by a Kelvin-Hclmholtz insta- 
bility, but also during the inspiral, most likely because the magnetic tension reduces 
the stellar tidal deformation for extremely large initial magnetic fields, Bq > 10 17 G. 
We quantify the influence of magnetic fields by computing the overlap, 0, between 
the waveforms produced during the inspiral by magnetized and unmagnetized bina- 
ries. We find that for any realistic magnetic field strength Bo < 10 14 G the overlap 
during the inspiral is O > 0.999 and is quite insensitive to the mass of the neutron 
stars. Only for unrealistically large magnetic fields like Bo ~ 10 17 G the overlap does 
decrease noticeably, becoming at our resolutions O < 0.76/0.67 for stars with baryon 
masses Mb ~ 1.4/1.6 M®, respectively. Because neutron stars are expected to merge 
with magnetic fields ~ 10 8 — 10 10 G and because present detectors are sensitive to 
O < 0.995, we conclude that it is very unlikely that the present detectors will be able 
to discern the presence of magnetic fields during the inspiral of neutron stars. 

Key words: relativity - gravitational waves - stars: neutron - binaries: general - 
magnetic fields - MHD 



1 INTRODUCTION 

Numerous astronomical observations suggest that large 
magnetic fields are associated with neutron stars (NSs). In- 
deed, evidence for the existence of binary NSs is obtained 
from binary pulsars, in which one or both NSs are seen to 
have a large magnetic field. In General Relativity such bi- 
nary systems cannot be stationary because they emit gravi- 
tational waves (GWs) which extract energy and angular mo- 
mentum from the binary, inducing it to inspiral and merge. 
During the final stages of the inspiral the GW emission is 
expected to be strong enough to be relevant for the detectors 
now operative at design sensitivities and it promises to pro- 
vide important informatio n on the equation of state (EOS) 
regulating the NS matter ijRead et al.ll2009h . In addition to 
their importance as sources of GWs, however, the merger 
of binary NSs is likely to provide important information on 
the physics of short gamma-ray bursts (GRBs). The coales- 
cence of the two NSs, in fact, gives rise, either promptly or 
after some interval, to a system composed of a torus orbit- 
ing a r ound a rapidly rotating black hole (BH) (jBaiotti et al.l 
|200S| ; lYamamoto et all 120081 ) . The complex plasma physics 



accompanying t his event is probably behin d the "engine" 
powering GRBs (|PirarJl2004 lMeszarodl2006h . 

There is little doubt, therefore, about the impor- 
tance of assessing the role played by magnetic fields in 
the inspiral and merger of binary NSs. Yet, determin- 
ing this accurately is a remarkably difficult task requiring 
the solution of the Einstein equations together with those 
of general-relativistic magnetohydrodynamics (GRMHD). 
So far, only two GRMHD simulations h ave been re- 
ported I Anderson et al.ll2008l ; iLiu et alj|2008l ). reaching dif- 
ferent conclusions about the importance of very strong mag- 
netic fields (B ~ 10 16 - 10 17 G). The aim of this Letter is 
to go beyond these qualitative estimates and provide a first 
quantitative measurement of the influence of magnetic fields 
on both the inspiral and the merger of magnetized NSs. 

By considering a large range of magnetic fields, which 
includes values more realistic than those used in the works 
cited above, and two different masses, we find that mag- 
netic fields generally grow after the merger, when the tur- 
bulent motions triggered during the merger by the Kelvin- 
Hclmholtz (KH) instability amplify any initial poloidal mag- 
netic field producing a toroidal one whose strength rapidly 
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becomes comparable to the poloidal one. In addition, we 
find that magnetic fields can, at least in principle, play a 
role already during the inspiral if sufficiently strong. This 
is most likely due to the magnetic tension, which decreases 
the NS deformability, increases the compactness, and thus 
delays the time of merger. In practice, however, the influ- 
ence of magnetic fields during the inspir al appears only for 
values ~ 10 1 7 G w hich are unrealistic (|Urpin et al.l 1 19981 ; 
lAbdolrahimil 120091 ). As a result, it is very unlikely that 
present detectors will be able to measure the presence of 
magnetic fields during the inspiral of NSs. Finally, we show 
that high-order numerical schemes are essential to draw ro- 
bust conclusions while lower-order methods incorrectly sug- 
gest that even strong magnetic fields have no influence at 
all. 



2 MATHEMATICAL AND NUMERICAL 
SETUP 

All the results presented here were computed by solv- 
ing the GRMHD equations in the ideal MHD approx- 
imation (i.e. assuming an infinite electrical conductiv- 
ity) and in dynamical spacetimes. The evolution of the 
spacetime was obtained using the CCATIE code, a three- 
dimensional finite-differencing code providing a solution 
of a conf ormally traceless fo r mulation of the E instein 
equations (jKoppitz et all 120071 ; IPollnev et all I2007I ). The 
GRMHD equations were instead solved usin g the Whisky 
code presented in iGiacomazzo fc Rezzollal (|2007l ). thus 
adopting a flux-conservative formulation of the GRMHD 
equations (|Ant6n et al.l 120061 ') and high-resolution shock- 
capturing schemes. In particular, we have computed the 
fluxes using the Harten-Lax -van Leer-Einfeldt (HLLE) ap- 
proximate Riemann solver (|Harten et al.l Il983h , while the 
reconstruction was made using the 3rd-order piecewise 
parabolic method (|Colella and Woodward! Il984l ). Further- 
more, to guarantee the divergence-free character of the 
MHD equations we have employ ed the flux-constrained- 
transport approach (|Totbl l2000i ). The code has been 
vali dated against a series of t ests in special relativ- 
ity llGiacomazzo fc Rezzollal 120061') and in full general rel- 
ativitv [see IGiacomazzo fc Rezzollal (|2007l )]. 

The system of GRMH D equations is closed by an EOS 
and, as discussed in detail in lBaiotti et alJ l)2008l ). the choice 
of the EOS plays a fundamental role in the post-merger dy- 
namics and significantly influences the survival time against 
gravitational collapse of the hyper-massive neutron star 
(HMNS) likely produced by the merger. It is hence impor- 
tant that special attention is paid to use EOSs that are 
physically realistic, as done in lOechslin and Jan 3 (120071 ) 
within a conformally-flat description of the fields and a sim- 
plified treatment of the hydrodynamics. Because we are here 
mostly concerned with computing a first quantitative esti- 
mate of the role played by magnetic fields rather than with a 
realistic description of the NS matter, we have employed the 
commonly used "ideal- fluid" EOS in which the pressure P 
is expressed as P = pe(F — 1), where p is the rest-mass den- 
sity, e is the specific internal energy and F is the adiabatic 
exponent. While simple, such EOS provides a reasonable 
approximation and we expect that the use of realistic EOSs 
would not change the main results of this work. 



All equations are solved on a Cartesian grid us- 
ing the vertex-centered mesh-refinement sche me provided 
by t he Carpet dr i ver llSchnetter et alJ 120041 ). Differently 
from iBaiotti et al.l (|2008t ) . we use here larger fixed refined 
grids rather than smaller moving ones. While computation- 
ally more expensive, this choice reduces the violations in 
the divergence of the magnetic field due to interpolations 
in the buffer zones between refinement levels. In this way, 
the divergence of the magnetic field on the finest grid (not 
including the buffer zones) is zero to machine precision. 

We have used five refinement levels with a 180-degree 
rotational symmetry around the z axis and a reflection sym- 
metry across the z — plane (in practice, we simulate only 
the region {x ^ 0, z ^ 0}). The finest grid has a resolution 
of h = 354.4 m and extends up to r = 44 km; the coarsest 
grid has h = 5.6704 km and extends up to r — 380 km. Our 
finest grid therefore contains both N Ss at all times and each 
NS is covered with « 80 :! points. In lAnderson et alJ (|2008h 
the finest grid had h = 0.46 km, th us with ~ 70 3 points 
across each star, and lLiu et al.l (|2008l ) had > 40 3 points. So 
our resolution is highe r than that in the ab ove works, but it 
is only half of that in lBaiotti et al.l l)2008l ). and it is barely 
sufficient to reach convergent results for the inspiral. 



3 INITIAL DATA. 



The i nitial data are the same as used in IBaiotti et al 



( | 200Sh and were produced by iTaniguchi fc Gourgoulhonl 



(|2002l ) with the multi-domain spectral-method code LORENE 
(http://www.lorene.obspm.fr). The initial solutions for 
the binaries are obtained assuming a quasi-circular orbit, 
an irrotational velocity field, and a conformally-flat spa- 
tial metric. The matter is modelled using a polytropic EOS 
P — Kp r with K = 123.6 and F = 2. Since no self-consistent 
solution is available yet for magnetized binaries, a poloidal 
magnetic field is added a-posteriori using the vector poten- 
tial A$ = vj 2 A b max (P — P cu t, 0)" s , where zu = *J x 2 + y 2 , 
Ab > parameterizes the strength of the magnetic field, P cu t 
defines where in the NS the magnetic field goes to zero, and 
n s determines the smoothness of the potential. The compo- 
nents of the magnetic field are then computed by taking the 
curl of the Cartesian components of A$ to enforce that the 
divergence of the magnetic field is zero to machine precision. 
Here we set P cut = 0.04max(P), and n s = 2 to enforce that 
both the ma gnetic field and i ts firs t derivative are zero at 
P — Pcut- In lAnderson et all (|2008f ) the magnetic field was 
built with an equivalent expression but with Pcut se t to the 
pressure in the atmosphere, and in iLiu et al.l <|2008h the ex- 
pression used is only slightly different, but P cut is set to be 
4% — 0.1% of max(P). Both above works set n a = 1. Note 
that the magnetic fields are confined at all times inside the 
NS matter and hence they cannot "repel" each ot her during 
the inspiral, as claimed in lAnderson et all (|2008l ) . 

Table Q] lists some of the properties of the eight equal- 
mass binaries studied here. In more detail, we consider two 
classes of binaries differing in the initial masses, i. e. binaries 
M1.45-B* (or low-mass), and binaries M1.62-B* (or high- 
mass). For each of these classes we take four different mag- 
netizations (indicated by the asterisk) so that, for instance, 
M1.45-B12 is a low-mass binary with a maximum initial 
magnetic field Bo = 1-97 x 10 12 G [the binaries with zero 
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Table 1. Properties of the eight equal-mass binaries considered. 
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(gcm~ 3 ) 




M1.45-B* 


or 1.97 x 10* 


14.3 


1.445 


2.681 


6.5083 x 10 49 


1.78 


15.2 


0.899 


4.58 X 10 14 


0.12 


M1.62-B* 


or 1.97 x 10* 


13.2 


1.625 


2.982 


7.7805 x 10 49 


1.85 


13.7 


0.931 


5.91 x 10 14 


0.14 



Notes. The different columns refer, respectively, to model name; maximum initial magnetic field Bo, where * is either (in which case 
-Bq = 0), 12,14 or 17; proper separation between the stellar centres d/M ADM ; baryon mass Mj, of each star; total ADM mass Af ADM ; 
angular momentum J; initial orbital angular velocity Co; mean coordinate radius r e along the line connecting the two stars; ratio of the 
polar to the equatorial coordinate radii r p /r e ; maximum rest-mass density p max ; compactness of the stars (M, / R)oo ■ 




t-r (ms) 



Figure 1. Top panel: £ = 2, m = 2 component of the /i+ polar- 
ization from binaries with different masses (thick line: 1.62 Mq; 
thin line: 1.45 Mq) and zero magnetic field. Bottom panel: the 
same but for binaries with an initial magnetic field Bq ~ 10 17 G. 

magne tic fields are the same as those evolved in lBaiotti et al.l 
(|2008t )]. In summary, we consider 8 different binaries that 
are either unmagnetized or with magnetic fields as large as 
Bo — 10 17 G. We reinforce a remark already made in the 
Introduction: all astronomical observations and theoretical 
considerations suggest that the magnetic fields in binary 
NS systems just before the merger are much smaller than 
~ 10 17 G. However, we do consider such l arge fields here, 
firstly to com pare with t hose studied by lAnderson et al.l 
l|2008h and by [[fiT et al.l (|2008l ). and secondly because by 
doing so we can determine important lower limits on the 
detectability of magnetic fields during the inspiral. 



4 GRAVITATIONAL WAVES AND OVERLAPS 

We postpone the discussion on the matter dynamics to a 
subsequent paper and concentrate here on the GW emission. 
A representative summary is offered in Fig. [T] which reports 
the I — 2, m = 2 component of the h+ polarization (mod- 
ulo a phase difference, h x shows the same behaviour). More 
specifically, Fig.[T]highlights the differences in the GWs from 
binaries with different masses, namely 1.62 M@ (thick line) 
and 1.45 Mq (thin line), when the initial magnetic field is ei- 
ther zero (top panel) or as high as ~ 10 17 G (bottom panel). 
Figure [1] clearly shows that when the NSs are not mag- 
netized, the high-mass binary has a larger-amplitude GW 
emission, it experiences an earlier merger and the HMNS 



at these resolutions collapses to a rapidly rotating BH after 
only ~ 4 ms, while the HMNS from the low-mass binary does 
not collapse. In contrast, when the NSs are initially magne- 
tized, the strong magnetic tension most likely reduces the 
tidal deformations and results in a delayed merger time (de- 
fined as the time when the maximum rest-ma ss density has 
a first significant minimum; cf. Fig 2 or 8 in iBaiotti et al.l 
Furthermore, the additional pressure support com- 
ing from the intense magnetic fields is such that neither the 
high-mass nor the low-mass binary collapse promptly to a 
BH over the ~ 15 ms of the simulations (see the bottom 
panel of Fig. for a comparison of GWs for high-mass bina- 
ries with different magnetic fields). Overall, Fig.[T]shows that 
magnetic fields have a strong impact on the GWs emitted 
after the merger but also during the inspiral, if sufficiently 
strong. If, on the other hand, the magnetic fields are more 
realistic, e.g. ~ 10 12 G, then differences appear only after 
the merger. For compactness, comparisons of this type can- 
not be presented here but will appear in a longer companion 
paper (Baiotti et al., in preparation). 

While generic, this behaviour depends sensitively on the 
strength of the initial magnetic field and there exists a crit- 
ical magnetic field below which the MHD effects during the 
inspiral are not important. In order to quantify this we have 
computed the overlap between two waveforms /i B1 , /i B2 from 
binaries with initial magnetic fields Bl, B2 as 

o[h B1 ,h m ]= - < h BiK 3 ) ^ =t (1) 

where {h B1 \h B2 ) is the scalar product and is defined as 

(h B1 \h B2 ) = ^ j~ ' df hB1 f*»> U) , (2) 

and h(f) is the Fourier transform of the GW h(t) and Sh(f) 
is the noise power spectral density of the detector (we have 
here considered LIGO). Clearly, waveforms that are very 
similar have O ~ 1. A general view is shown in Fig. [2] which 
reports the overlaps between the unmagnetized binaries and 
binaries with different magnetizations, i.e. O[h B0 , h g ], for 
the two masses considered here (top and bottom panels, re- 
spectively). Note that the overlap is relative to the inspiral 
only [i. e. the integral @ is cut off at the orbital frequency 
at merger] since this is the phase for which our results are 
convergent (becoming only consistent after the merger as a 
result of the development of turbulence) and because the 
post-merger evolution can only further decrease O. It is evi- 
dent that for the high-mass binary (top panel) the influence 
of the magnetic field is noticeable only for very large mag- 
netic fields (O ~ 0.999 for B ~ 10 14 G and O ~ 0.668 
for Bq ~ 10 17 G) . This is true also for the low-mass binary 
(bottom panel) whose smaller compactness, however, leads 
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Figure 2. Overlap in inspiral waveforms from binaries with dif- 
ferent magnetization; the top and bottom panels refer to the high- 
and low-mass binaries, respectively. 

to larger overlaps (i.e. O ~ 0.761 for Bo — 10 17 G). In view 
of these results, of the fact that NSs just prior t o merger are 
expected to have magnet ic fields ~ 10 8 — 10 10 G |Urpin et al.l 
1 19981 : I Abdolrahimill2009ri , and that presen t detectors are sen- 
sitive to O < 0.995 (|Lindblom et al-boOST ), we conclude that 
it is very unlikely that the presence of realistic magnetic 
fields can be detected during the inspiral. 



5 MAGNETIC FIELD AMPLIFICATION 

As discussed in detail in iBaiotti et al.l |2008l ). during the 
merger a shear layer develops in the region where the two 
NSs enter in contact. Across this layer, the tangential com- 
ponents of the velocity are discontinuous and this leads to 
the development of a KH insta bility and thus to the produc- 
tion of vortices [cf. Fig. 16 of IBaiotti et al.l (|200ot )]. When 
poloidal magnetic fields are present, this hydrodynamical 
instability can lead to exponentially growing toroidal mag- 
netic fields, thus increasing the energy stored in magnetic 
fields. Th is mechanism, already obs erved in Newtonian sim- 
ulations l)Price and Rosswod l2006t ) but not before in gen- 
eral relativistic ones, is likely to be important for explain- 
ing the physics powering short GRBs. Taking M1.62-B12 
as a reference, Fig. [3] shows the evolution of the maxima 
of the magnetic field |J3| = (B i B i ) 1/2 (thick line), and of 
its toroidal \B T \ (dashed line), and poloidal |B P | (thin line) 
components. A vertical dotted line marks the merger, occur- 
ring « 1 ms after the KH instability has started developing. 
Clearly, as long as the KH instability is active, the toroidal 
magnetic field is amplified exponentially, until it reaches val- 
ues comparable to the poloidal o ne (this is different and mor e 
reasonable than what found by IPrice arid" Rosswod (|2006h , 
where the magnetic field reached energy equipartition val- 
ues). Note that the magnetic field grows considerably also 
when the HMNS collapses to a BH as a result of magnetic- 
flux conservation in the collapsing NS matter . 

As discuss ed in IPrice and Rosswod (|2006l ) and 

in IBaiotti et al.l (|200&f ). much in the development of the 



Figure 3. Evolution of the maxima of the magnetic field \B\ = 
(B l Bi) 1/l2 (thick solid line) and of its toroidal \B T \ (dashed line) 
and poloidal \B P \ (thin solid line) components for the high-mass 
case with So — 10 12 G. The dotted line marks the merger time. 

KH instability and in the subsequent magnetic field ampli- 
fication depends on the resolution used. A detailed study 
of the turbulent regime and magnetic field amplification 
produced by the merger is extremely challenging and re- 
quires resolutions well above the ones that can be afforded 
now in GRMHD simulations. Nevertheless, on the basis 
of preliminary investigations with different resolutions, we 
expect the behaviour in Fig. [3] to be qualitatively correct 
and, hence, that as long as the KH is active, the poloidal 
magnetic field is coiled into a toroidal one, increasing it 
exponentially to values comparable with the poloidal one. 
When equipartition among the components is reached, 
the large magnetic tension suppresses the KH instability, 
preventing a further growth of the toroidal magnetic field. 
An analysis of this process will be presented elsewhere. 



6 THE IMPORTANCE OF HIGH-ORDER 
METHODS 

As shown in IBaiotti et al.l < |2008h . the use of reconstruction 
schemes of sufficiently high order is essential for a correct 
calculation of the GW signal. This is stressed also in Fig. [4] 
which presents a comparison in the GW emission from the 
high-mass binary for evolutions made using either a 2nd- 
order MINMOD scheme (top panel) or a 3rd-order PPM 
one (bottom panel) with the same grid structure and res- 
olution used in the previous runs. In both cases a thick 
line refers to the unmagnetized binary while a thin line to 
the binary with Bo ~ 10 17 G. Clearly, while the evolutions 
with MINMOD show only minimal differences between the 
magnetized and unmagnetized case (O = 0.9994 over the 
whole waveform), the evolutions using PPM show consid- 
erable differences (O = 0.6500), both during the inspiral 
and after the merger. More precisely, although we use ex- 
actly the same initial data, the binaries evolved with MIN- 
MOD merge almost two orbits earlier than those evolved 
with PPM [cf. vertical dotted lines in Fig. H]. Additionally, 
the unmagnetized binary evolved with MINMOD does not 
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Figure 4. GW emission from the the high-mass binary evolved 
using either a 2nd-order scheme (top panel) or a 3rd-order one 
(bottom panel). A thick (thin) line refers to a binary with Bo = 
0G (10 17 G), while the vertical lines indicate the time of merger. 

collapse to a BH, in contrast to what happens when using 
PPM at these resolutions. These differences are due to the 
numerical dissipation of the 2nd-order method, which is in- 
adequate at th ese resolu ti ons. This could explain why the 
calculations in iLiu et al.l (|2008l ), where a 2nd-order recon- 
struction and a lower resolution were used, show only small 
differences between unmagnetized and magnetized binaries. 



7 CONCLUSIONS 

We have presented accurate simulations of the inspiral and 
merger of magnetized NSs and found that magnetic fields 
have an impact both during the inspiral and after the merger 
but only if sufficiently strong. Comparing waveforms for dif- 
ferent magnetizations we have found that for Bo < 10 14 G, 
the overlap O > 0.999 and is quite insensitive to the mass 
of the NSs. Only for unrealistically large magnetic fields 
(Bo — 10 17 G), the overlap decreases noticeably, becom- 
ing at our resolutions O < 0.76/0.67 for stars with baryon 
masses M& ~ 1.4/1.6 M©, respectively. Since the magnetic 
fields in NSs just prior to merger are expected to be rather 
small (~ 10 s - 10 10 G) (jUrpin et all Il998l ; lAbdolrahiml 
l2009h . we conclude that it is very unlikely that the present 
detectors will be able to measure the presence of magnetic 
fields during the inspiral. Magnetic fields could be however 
detectable after the merger and hence in th e part of the 
spectr um at frequencies > 2kHz [cf. spectra in lBaiotti et al.1 

Another important result discussed here is the evidence 
that a KH instability develops during the merger, leading 
to the exponential growth of a toroidal magnetic field, the 
strength of which becomes comparable with the poloidal 
one. This additional magnetic field can modify the structure 
of the HMNS and decrease the overlap after the merger. Fi- 
nally, we have provided concrete evidence that high-order 
methods are essential to draw robust conclusions and that 
instead lower-order methods incorrectly suggest that mag- 
netic fields have no influence at all. 
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